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Abstract 

We introduce a solvable quantum antiferromagnetic model. The model, with Ising spins in 
a transverse field, has infinite range antiferromagnetic interactions and random fields on each 
site following an arbitrary distribution. As is well-known, frustration in the random field Ising 
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frustration and we report some analytical results for it. 
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I. INTRODUCTION 



With the realization in the mid last century that the Neel state cannot be the ground 
state (not even an eigenstate) of a quantum Heisenberg antiferromagnet, considerable effort 
has gone in search of, and in understanding the nature of, the ground state of such and 

n 

similar quantum antiferromagnet Since early 1960s, quantum spin systems described 
by Ising model in a transverse tunneling field was investigated extensively; particularly be- 
cause of easy mapping of the quantum system to its equivalent classical system and some 
cases of exact solubility However there have, so far, been very few soluble models with 
antiferromagnetic interactions. It is well-known that the Ising model with long range inter- 
actions is solved exactly, even if the system has some special kind of quenched disorder, like 
in Sherrington-Kirkpatrick model of spin glasses. The number of degenerate states there 
can be estimated to be 0(2 N / 2 ), which is larger than that of the above mentioned model 
(0(2 a28743Ar )). However, it is not so easy to consider the antiferromagnetic version of the 
model due to a lack of sub-lattice to capture the Neel ordering at low temperatures. In 
this paper, we introduce and study a solvable quantum antiferromagnetic model. In our 
model system each spin is influenced by the infinite range antiferromagnetic interactions 
in a transverse field. We also consider the case under the Gaussian or the binary random 
fields. By introducing two sub-groups of the spin system, we describe the system by means 
of the effective single spin Hamiltonian which is derived by the Trotter decomposition 3 
and Hubberd-Stratonovich transformation j^, and we solve the model analytically. A pre- 
liminary study, for the case without random fields, was reported earlier jsj]. 

This paper is organized as follows. In the next section, we introduce our model system 
and write down the general formula of the averaged free energy density. In section 3, to 
check the validity of our analysis, we compare our result with the previous well-known result 
which was obtained by mean- field approximations [6]. In section 4, we consider the system 
under the Gaussian and the binary random fields on site and derive the equations of states 
and evaluate them. We then obtain the phase diagrams. Section 5 gives a summary. 
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II. THE MODEL SYSTEM AND ITS ANALYSIS 



In order to capture the Neel ordering below the critical temperature Tjy or the am- 
plitude of transverse field T N , we divide spins S into the sub group A : S^ A ' = 
(S X ' Z ' {A \ S* MA \ ■■■ , S X /' (A) ) and the sub group B : S {B) = (S*' z ' iB \ S X ' Z ' (B) , ■■■ , S X /' {B) ), 
which are corresponding to virtual sub-lattice A and B. Then the system is described by the 
following effective Hamiltonian 

H <?) = jj E s * ,iA)sZ / B) - E r <E *r (0 - h E E k# (0 > w 

ij Z=A,B i 2=A,S i 

where S*'^ A ' B ', Sf'^ A ' B ^ are x and z components of the Pauli matrix : 

and ft. = (/ii, /12, • • • , /ijv) is a vector of the random fields on site and h means the strength 
of the random field. Ta and Tb are amplitudes of the transverse fields in the sub groups A 
and B. The main advantage of this Hamiltonian is that it can be recast exactly to that of a 
single spin in an effective field. 

Using the Suzuki- Trotter formalism one can express the quenched-variable h- 
dependent partition function as 

Z(h) = lim tro <? e"^) (3) 



with 

P H W = OWU E 1 E ^ + E S ik ] ( ~~ oWM E 1 E ^ ~ E ^ 



where 7a,b = |logcoth \^-^- \ ■ By using the Hubberd-Stratonovich transformation 
the field ft-dependent partition function is written by means of the saddle point technique, 
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in the limit N — » oo, as 
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(5) 



i=l 



where H^hj) is the effective single spin Hamiltonian and is given by 



M 



(6) 

Here we should keep in mind that the integrals with respect to m k + and m k _ are evaluated 
at the saddle points in the limit of N — > oo, namely 



m 



4 E C } + ^E = -( M " z + 



+ 



iV 



iV 



to. 



(7) 
(8) 



i i 

where to+ and m/L are related to the magnetizations of z-component for group A and B, 
namely M k / and M^ z . Thus the free energy F(h) = -f3^ log Z(h) of the system is now 
written by 

N 

F ^ = E«) 2 + E( m -) 2 " loglltr^ ^exp [-^{K)] . (9) 

k k i=l 

Taking into account the symmetry of the system, to+ = m + ,m k _ = to_ for all k, i.e., 
the so-called static approximation holds good naturally; the fluctuations due to, say, two- 
spin correlations (including entanglements) vanishes in thermodynamic limit. Thus the 
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/i-dependent free energy leads to 

N 

F(h) = ~{M% + Ml) 2 + *±{M\ - M|) 2 - (3-Ho g \{Z^Z^ (10) 

i=i 

with 

i 

where we used the relations (JJJ) and (JHJ): m + +m_ = — 2M#, m + — m_ = —2M Z A . It should be 
noted that the above free energy still depends on the fields h. To obtain the /i-independent 
averaged free energy F, we should evaluate the following quantity 



/oo 
P(h)F(h)dh, (12) 
-oo 



where P(h) = P(hi,--- ,hw) is a joint distribution of the random fields and we defined 
dh = dhi ■ ■ ■ dhpj. If we assume that the random field hi for each site % is uncorrelated (not 
influenced by other h/s; j ^ i) then 

P(h) = P{h u ■ ■ ■ , h N ) = P(hi) ■ ■ ■ P(h N ) = 11 P{h t ), (13) 

i 

there in the average in (|12j) . giving the free energy (per spin) 

fi = -JMfM^ x -^ I" dhP(h) log 2 cosh (3^/{2 JM* - hhy + T 2 +1 . (14) 

l=A,B 

Here we used the fact that the 2x2 matrix H (with the elements {H)n — —(H) 2 2 = 
a, (H)i2 = (H)2i = b appearing in the of exponent of equation (JHJ)) has eigen values 



±\/a 2 + b 2 . We also should keep in mind that in the sum with respect to I, A + 1 = 
B,B + 1 = A; hence \A — B\ = 1 is satisfied. Hereafter, we use this relation for the sum 
with respect to the label /. The magnetizations for two sub-lattices M\ and M% now obey 
the following saddle point equations 

Mi= ( dhPCh) — ^ 2JMl+1 = kk )=t&nhpJ(2JM l z +1 - hh) 2 + Tf, (I = A, B) (15) 
-l(2JMf +l -hh) 2 + T 2 



where it should be noted that A — B = 1 and the appropriate expression for magnetization 
of each sub-lattice can be derived by equating df /dM^ and df jdM z A to zero respectively. 
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III. ANALYSIS UNDER UNIFORM FIELD 



We first consider the case of uniform field P|, i.e., in (JT%j) or 1)150 . /i — > 1. In this limit, 
by using the fact dhP(h) = 1, the free energy density / reads 



/ = E ^ = ~ JM ? M W - /?- 1 log2cosh/?^(2JMf - hf + rf +1 (16) 
and the saddle point equations with respect to M A and M B are obtained as follows. 



(—2TM Z 4- h) i 

Mf= 1 m J =tanh/3,/(-2JMf +1 + /i) 2 + r ; 2 , {l = A,B). (17) 

;-2 JMf +1 + hy + r 2 v 



Before we investigate the quantum effects, we check the classical case, that is Ta = ^ b = 0. 
Then, the above saddle point equations are reduced to 

Mf = tanh/3(-2 JMf +l + h), (I = A, B). (18) 

In order to determine the Neel temperature, we expand the equations around Ma, Mb — 
for h = and obtain M Z A ~ —2J(3M B ,M B ~ —2J(3M\. From these linearized equations, 
we find that only possible solution for the case of M Z A = M B = M z is M z = 0; whereas with 
Neel ordering a finite value of M z can be obtain: M Z A = —M B = —M z . This gives the Neel 
temperature T/v = P^ 1 = 2J. The linear susceptibilities xa and xb are then evaluated as 

A ^-oo dh cosh 2 f3( -2 JM[ +1 ) V ; V ' 

The behavior beyond the Neel temperature T > is determined by the condition M Z A ~ 
0, M Z B ~ 0, i.e., xa = - 2 J^b) and xb = - 2 J^a)- This leads to 

X = Tf^Tf- (20) 
l + In 

where we defined X = Xa + Xb- Therefore, in the limit of T — >• T/v = J, the linear 
susceptibility x converges to x ~* 1/2 J. On the other hand, below the Neel temperature 
T < Tn, we find the solution for the Neel ordering : M A = —M B = —M z ^ for h = 0. 
This condition should be satisfied for (cf. Eqn. (31)) 

x = 2 1 2 (21) 

2/3 J + cosh (3(2 J M z ) v ; 

and at the critical point T = Tjv, M = 0, the susceptibility takes x — 2/(2 J + T/v) = 1/2 J. 
In Fig. [T] (left) we plot the shape of the susceptibility x as a function of T. From this 




FIG. 1: The longitudinal magnetization M z as a function of T and T for uniform systems (the left panel). 
The critical points of the second order phase transition are given by TV = Tfj = 1/2 J. The right panel 
shows the corresponding susceptibility \ as a function of T and T. 

figure we find that the susceptibility has a cusp, instead of the divergence as observed in the 
ferromagnetic systems, at the critical temperature T = Tjv, as observed in the analysis for 
finite range models by using mean-field approximations (|. As mentioned earlier, this model 
being recastable exactly to a (classical) single-spin in an effective field (i.e., the mean-field 
approximation for the partition function being exact), there is no scope of any non-trivial 
wave-function with entangled neighbouring spins. 

We next consider the quantum case Ta,T b ^ 0. For simplicity, we consider the case 
of the symmetric transverse field, namely, = Tb = T. In order to consider the pure 
quantum effects, we take the limit of (3 — > oo and deal with the following coupled equations 

-2 JM, Z , , + h 

Mf = 1+1 (I = A,B). 22 

It is important for us to bear in mind that for J > 0, the above equations have a solution 
M z = if M A = M B and a solution M z ^ if M A = -M B , as expected. To determine 
the critical transverse field, we expand the above equations around Mj[, M B ~ for h = 
as M% ~ -2JM Z B /T and M Z B ~ -2JM Z A /T. This gives the critical point T N = 23. The 



spontaneous sub- lattice order M Z A or M B vanishes at the Neel phase boundary Tn(T). Deep 
inside the antiferromagnetic phase (at (3 — ► oo,T — ► 0, h = 0), = 1 = — M B and the 
free energy density / can be expressed as / = l//31og[l + exp(— /3A(T))], the specific heat 
d 2 f/dT 2 will have a variation like exp[— /3A(r)] like that of a two level system with a gap 



A(r) = \/4J 2 + T 2 here. This is the exact magnitude of the gap in the magnon spectrum 
of this long range transverse Ising antiferromagnet. 
The susceptibilities are given by 

dMf r 2 (i-2j Xm ) 
^ = i im o^r = [(2JM^ 1 ) 2 + r 2 ]3/ 2 ' ( / = A ^ < 23 > 

Then, the behavior beyond the critical amplitude of the transverse field is determined 
by the condition M|, M| ~ 0, namely, X a = (1 - 2Jxb)/T, Xfl = (1 - 2Jx j4 )/r. This leads 
to 

2 

X = Xa + Xb = p , -p • (24) 

1 + 1 N 

On the other hand, below the critical amplitude of the transverse field T^v, we set M Z A = 
M Z B = —M z and obtain 

x ~ [(2JM) 2 + r 2 ] 3 / 2 + 2jr 2 ' ( ^ 

At the critical point T N [M z = 0), the susceptibility is given as x — ^/(Xn + 2 J) = 1/2 J. 
Around T ~ 0, the susceptibility behaves as x = 2r 2 /(2J) 3 . In Fig. U (right), we plot x as 
a function of T (and also T). From this figure, we find that the susceptibility has a cusp at 
the critical amplitude of the tunneling field. 

We next investigate the transverse component of the susceptibility. The magnetization of 
the transverse direction are calculated the derivative of the free energy density with respect 
to the amplitudes of the transverse field T A , T B . 

M? = ^ = — Tl =tanh/3J(2JMr +1 + fr) 2 + r 2 , (/ = A, B) (26) 

dTl y/(-2 JM? +1 + hf + r 2 V 

At the ground state, these coupled equations are simplified as follows. 

Mf = , F * {I = A, B) (27) 

-2Mfa + hf + r 2 



s 



In para-magnetic phase is specified by M A = M B = and this gives M A = M B 



1. 



On the other hand, antiferromagnetic phase for h — 0, we obtain from 
~JM Z B / V /(2JM Z B ) 2 + T 2 AJ M Z B = -JM z A /^{2JM z A y + T 2 B . Hence, 



as M Z A 



Mf 



(28) 



for the antiferromagnetic case M A = —M B = —M. Therefore, the susceptibilities of the 
transverse direction lead to 



X" 



for Ta, r# > Tn and 



Xi 



X 



dT, 



0, (l = A,B) 



(29) 



Xa = Xb 



J 



(30) 



M A = M B for M z 



for Ta,^b < Tat- We plot the transverse magnetization M x 
M B in Fig. |21 We next consider the quantum antiferromagnetic system under random fields. 




FIG. 2: Longitudinal magnetization M x and susceptibility \ A 



IV. ANALYSIS UNDER RANDOM FIELDS 



In the previous section, we investigated the critical phenomena for spatially uniform 
systems. It has been conjectured that fluctuation in the random (longitudinal) field Ising 
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model (RFIM) gives rise to a many valley structure in the configuration space, similar to the 
case in spin glasses. The study of the longitudinal random field transverse field Ising model 
with ferromagnetic uniform interactions has already been made jjj. There seems to be no 
reported analytic research for the RFIM with antiferromagnetic interactions in transverse 
field. It should be interesting to investigate the competition between two different kinds 
of frustration; frustration due to the quenched disorder and the frustration induced by 
the antiferromagnetic interactions. In this section, we investigate quantum systems under 
random on-site longitudinal fields. In this paper, we consider the following two cases of the 
random field distributions : 



P g (h) = f[-L-ex P -±^( hi -h ) 2 = f[P g (hi) (31) 

8=1 v L J 

N N 

and P b (/i) =Y[{08(h i -h o ) + (l-0)6(h i + h o )} = l[P b (h i ), (32) 

i=l i=l 

where the bias factor of the binary random field 9 takes < 9 < 1. For these distributions, 
the free energy densities f 9t b = dhP 9t b(h)F(h) become 



f 9 = E /' ( 33 ) 

l=A,B 

ff = --JMf.,Mf - 1 1 I Dx log 2 cosh pJ{ -2 JMf +1 + h(ax + h )} 2 + Ff, (34) 

J — oo 

for the Gaussian random field P g (h), where Dx = -J^e~ x ' 2 ^ 2 dx and 

h = E fi (35) 



=A,B 



= -JMf +1 Mf - p- x 6 log 2 cosh (3 J {-2Mf +1 J + hh } 2 + 



- (3- 1 {l-9)\og2cosh(3^{2M l z +1 J + hh } 2 + Tf, (36) 

for the binary random field Pb{h). In following, we investigate the critical phenomena given 
by these free energy densities for these two cases. 
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A. The Gaussian random field 



For the Gaussian random field (|3~T|) the saddle point equations are given by the derivative 
of the free energy density fa with respect to M\ and M B . Then we have 

f°° \-2JM z , + Max + ho)} n ; ^ o 

Mf = / Dx 1 1+1 v Jf tanh (3J{ -2 JMf +1 + /i(<xx + /i )} 2 + r 2 

^/{-2JMf +1 + /i( ( xx + /io)} 2 + r2 V 

(37) 

for I = A,B. At the ground state {(3 — > oo) these equations of states are simplified as follows 



{-2JMf +1 + %z + /i )} 2 + r 



For the possible choice M Z A = M z = —M^ to detect the Neel ordering, we solve the above 
equations for the symmetric transverse fields = Tb = T numerically. In Fig. El (left), we 
plot the case of the center of the Gaussian, h takes h Q = and the deviation of the Gaussian 
a = 0.5 and 1.5. From this figure, we find that the system undergoes second-order phase 
transition at the critical amplitude of the transverse field from the behavior of M z . If the 
phase transition is first order for the case of the center ho of the Gaussian (jUlj) is zero, we 
can expand the saddle point equation with respect to M Z A under the condition = Tb = T 
and M% = -M% = —M z as 

M z = C X M Z - C 3 (M 2 ) 3 + £>((M 2 ) 5 ) (39) 

with 

^i = 2JF 2 f — g^— , C 3 = AJ^ Dx ^°f + ZL (40) 



!_ O0 [(hcJxy + T 2 ] 3 / 2, 6 7_ 00 [(hax) 2 + r 2 ] 7 /2- 

The phase boundary T(a) of the continuous transition between the Neel and the paramag- 
netic phases is obtained for a given set of the parameters (J, ho, h) by the condition a = 1 
namely 

Dx N - 1/2 



2J loo[(^) 2 + r 2 ]3/ 2 J • (41) 

The second order phase transition is observed for C3 > 0, whereas a first order phase 
transition is found for C3 < 0. In Fig. El (right), we plot the boundaries between the Neel 
and the paramagnetic phases for both quantum and classical systems. In this plot, we set 
(J, h, ho) = (1, 1, 0). In the left panel of Fig. 01 we plot the factor C 3 of the third order of 
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the expansion of the magnetization M z as a function of a. In this plot we substituted the 
solution of the boundary (j4~Tj) for a given a into C3. From this panel we find, that from the 
value C3 = 4(J/r) 3 > at a = 0, that C3 decreases and takes its positive minima at just 
below the critical point a c , and beyond the critical point C3 increases again. We therefore 
conclude that the phase transition is always second order. We might see this from the 




r cr 



FIG. 3: The right panel shows phase boundaries between antifcrromagnetic and paramagnetic phases for 
classical and quantum systems for (J,h,ho) = (1, 1, 1). The left panel shows the variation of longitudinal 
and transverse magnetizations M z and M x with V. 

argument below. In the limit of F — > 0, the equation of state for M B = —M\ = —M z 
is simplified to 

M z = J Dx sgn (2 JM Z + hax) = 1 - 2H I— — J . (42) 

In Fig. 0] (right panel), we plot the solution of ()42j) for several values of h. We found that 
for cr — > 00, H{2JM z /ha) -> 1/2 and M z -> 0, whereas, for a -> 0, H{2JM z /ha) -> and 
AP — > 1. By expanding ()42j) with respect to M z up to the first order, we obtain the critical 
point a c as 
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FIG. 4: The longitudinal magnetization M z as a function of a for several values of ho for the case of T = 
(right panel). Each line of a solution of the equation 142f) . The left panel shows the cr-dependence of the 
factor C3 of the third order of the magnetization near the critical point. 

The magnetization varies continuously near this critical point (of the second order phase 
transition): M z = ^12a/ira c (l - a/a c ) 1 / 2 . 

We next evaluate the x-component (the transverse component) of the magnetizations M\ 
and Mg. These are calculated at the ground state (3 — > 00 as 

Mr = ^ = f , (i = A, B) (44, 



dT, 



{-2JM*+h{ax + h Q )} 2 + Y 



In Fig. El (left panel), we plot the results. It should be noted that in the limit of Y — ► 00, 

M x saturates M x ~ J^DxiT/T) = 1. 
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B. The binary random fields 



For the binary random fields, we obtain the saddle point equations by taking the derivative 
of the free energy density fs with respect to M\ and M| as follows: 

0(2JMf +1 - hh ) 



Mf 



{2JM X 



l+i 



hhn 



r? 



tanh/3J(2JM? +1 - hh ) 2 + rf 



9)(2JM z +1 + hh ) 



'(2JMf +1 + hh y + r 2 

At the ground state j3 — * oo, these equations are simplified to 



tanh (5 J (2 JM[ +1 + hh ) 2 + T 2 , (l = A,B). (45) 



Mf 



[l-6){2JMf +1 +hh ) 



9(2JM[ +1 - hh ) 



(2JMf +1 - hh ) 2 + T 2 J(2JM? +1 + hh ) 2 + T 



(l = A,B). (46) 



We solve the above equations numerically and plot it in Fig. El (left panel). From this figure, 
we find that the system undergoes first order phase transition when the value of ho is larger 
than same critical point h^. Whereas, for small value of ho < h^, the phase transition is the 
second order. 

In following, we determine the tri-critical point (h^,T c ). If the transition is continuous, 
we can expand the saddle point equation for M\ under the condition = T B = T and 



M Z B = -MX = —M z as follows: 



M z 



C + dM z + C 2 (M Z Y + C 3 (M Z Y + 0((M 



rz\4:\ 



where we defined 



Co 
C 3 



2(9 - l)hho 
^(hh ) 2 + T 2 
2JT 2 

{{hho) 2 + v 2 yi 2 

-2J 2 hh (29 - 1) 



Y 4 -AY 2 (hho 



r 2 + 4(hh ) 2 
_{(hho) 2 + T 2 } 5 / 2 

2 



(47) 

(48) 
(49) 
(50) 
(51) 



.{(/^ ) 2 + r 2 } 7 /2_ 

Therefore, if the distribution (|32~j) is symmetric (i.e., 9 = 1/2), the factors C and C 2 vanish 
and the magnetization behaves as 



M z = C X M Z + C 3 (M z f + 0((M z f). 



(52) 
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From this expression, we find that a second order phase transition is found when the condi- 
tion C\ = 1 and C3 < holds. On the other hand, a first order phase transition is observed 
for C 1 = 1 and C 3 > 0. Therefore, the point (/^,r c ) = (( J//i)(4/5) 3 / 2 , 2J(4/5) 3 / 2 ), which 
is determine by C\ — 1,6*3 — 0, corresponds to a tri-critical point on the phase boundary. 
In Fig. |5] (right panel) we plot these phase boundaries. We should notice that the critical 
point T c is independent of h. We find that for ho > h^, the transition from the symmetry 
breaking phase to symmetric phase is first order. To compare this result with the case of 
the Gaussian random field, we consider the limit r — > in (HBI). We obtain 



Mf = 9 sgn(2 JM? +1 + hh ) + (1 - 0) sgn(2 JMf +1 - hh ), (I = A,B). (53) 

To detect the transition point between the Neel and paramagnetic phases, we set Mf = 
M z = —Mf +l and 9 = 1/2 for simplicity. We then have 

2M Z = sgn(2 JM Z + hh ) + sgn(2JM z - hh ). (54) 

Apparently, M z takes values 1 or and the critical point of the first order phase transition 
is determined by 2 J — hh = 0, i.e., = 2J/h. This point is observed on the crossing 
point on the h -a.xis in Fig. (right panel). On the other hand, as we saw in Fig. 0] (left 
panel), the magnetization M z for the Gaussian random field drops gradually and the phase 
transition is second order even if there is no quantum fluctuation T = at the ground state 
(at (3 —>■ 00). This is a reason why the order-disorder phase transition in the Gaussian 
random field Ising model is always first order and it does not depend on T or a. Finally, we 
calculate the transverse magnetization M\ and Mg. From the derivative of the free energy 
density Jb with respect to and r#, we obtain 



= S = tanh/y{-2JMf +1 + /^ } 2 + r 2 

9r < M^JM^ + hhoY + Tl, 



!l 0)1 1 :tanh/5,/{2JMf +1 + hh } 2 + Tf, (l = A,B). (55) 



{2JM Z +1 + hh } 2 + r 2 
At the ground state, these equations are simplified as 

Mf= , S +- {l ~ e)Vl (l = A,B). (56) 



'{-2JM z +l + hh } 2 + T 2 J{2JMf +x + hh } 2 + r 2 

In Fig. El (left panel) we plot the transverse magnetization for case of the symmetric ampli- 
tude of the tunneling field = T B = T as a function of T. Obviously, for large T, we found 
M x = 1 from (jHU). 
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FIG. 5: The longitudinal and the transverse magnetizations, M z and M x as a function of T for the case of 
the binary random fields (right panel). The left panel shows phase boundaries between antiferromagnctic 
and paramagnetic phases. The dots represent tri-critical points (/iQ,r c ) = ((l//i)(4/5) 3/ ' 2 , 2(4/5) 3 / 2 ) we set 
J = 1. 

As may be noted, in the case where the random quenched fields are symmetrically dis- 
tributed about its zero value, the effective sub-lattice symmetry could be utilized to reduce 
the whole problem to that of a ferromagnet. 



V. SUMMARY 



In this paper, we proposed an analytically solvable quantum antiferromagnetic Ising 
model. Because of Ising anisotropy and long-range interactions, it has solvable Neel-like 
ground and other stateproperties. In view of the extensive recent studies in quantum 
a„tif_net S (Q-rf), particular* in q uantu m lsin g a„tif_net S (Q.fl), this 
kind of analysis should be of considerable importance. 

In the analysis of spatially uniform system, we found the Neel order below the tunneling 
field Ttv and show that the linear susceptibility has a cusp variation around that critical T^. 
It may be mentioned that a similar behavior in the half-filled Hubberd model was observed 



16 



earlier 




For this uniform spin system, the free energy density (JTfij) gives the ground state energy 
in the zero temperature limit and it also gives the low temperature behavior of the specific 
heat, the exponential variation of which gives the precise gap mag nitude A(= V4J 2 + T 2 ) 
in the excitation spectrum of the system. It may be noted that, because of the restricted 
(Ising) symmetry and the infinite dimensionality (long range interaction) involved, there 
need not be any conflict with the Haldane conjecture. Although our entire analysis has been 
for spin- 1/2 (Ising) case, because of the reduction of the effective Hamiltonian © to that of 
a single spin in an effective vector field, the results can be easily generalized for higher values 
of the spin S. No qualitative change is observed. The order-disorder transition in the model 
can be driven both by thermal fluctuations (increasing T) or by the quantum fluctuations 
(increasing T). This transition in the model has been investigated studying the behaviors 
of the (random sub-lattice) staggered magnetization and the (longitudinal and transverse) 
susceptibilities. No quantum phase transition, where the gap A vanishes, is observed in the 
model, unlike in the one dimensional transverse Ising antiferromagnets. 

By analysis of the disordered system as the random field Ising model in a transverse 
field, we found that the order of the phase transition changes at a tri-critical point. These 
conditions are obtained analytically for both the Gaussian random fields and the binary 
random fields. We believe, analysis of such model systems might provide some insights also 
for the quantum antiferromagnets with short range interactions. 
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